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' Abstract. We calculate light-cone distribution amplitudes for non-relativistic bound states, includ- 

ing radiative corrections from relativistic gluon exchange to first order in the strong coupling con- 
stant. Our results apply to hard exclusive reactions with non-relativistic bound states in the QCD 
factorization approach like, for instance. Be r\Av or e+e^ J /^r\c. They also serve as a toy 
model for light-cone distribution amplitudes of light mesons or heavy B and D mesons. 
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INTRODUCTION 



Exclusive hadron reactions with large momentum transfer involve strong interaction 
dynamics at very different momentum scales. In cases where the hard-scattering process 
is dominated by light-like distances, the long-distance hadronic information is given in 
terms of so-called light-cone distribution amplitudes (LCDAs) which are defined from 
\ hadron-to-vacuum matrix elements of non-local operators with quark and gluon field 

O ■ operators separated along the light-cone. Representing universal hadronic properties, 

LCDAs can either be extracted from experimental data, or they have to be constrained by 
^ ! non-perturbative methods. The most extensively studied and probably best understood 

^ I case is the leading-twist pion LCDA, for which experimental constraints from the n — y 

transition form factor [1], as well as estimates for the lowest moments from QCD sum 
rules [0, Hi |4J and lattice QCD [5J exist. On the other hand, our knowledge on LCDAs 
for heavy B mesons [SItIjIII], and even more so for heavy quarkonia [@, [I^Ij had been 
relatively poor until recently. 

The situation becomes somewhat simpler, if the hadron under consideration can be 
approximated as a non-relativistic bound state of two sufficiently heavy quarks. In this 
case we expect exclusive matrix elements - like transition form factors iflTIl and, in 
particular, the LCDAs - to be calculable perturbatively, since the quark masses provide 
an intrinsic physical infrared regulator. In these proceedings we report about results 



from [tl2D, where we have calculated the LCDAs for non-relativistic meson bound states 
including relativistic QCD corrections to first order in the strong coupling constant at 
the non-relativistic matching scale which is set by the mass of the lighter quark in the 
hadron. 



LIGHT-CONE DISTRIBUTION AMPLITUDES 



The wave function for a non-relativistic (NR) bound state of a quark and an antiquark 
can be obtained from the solution of the Schrodinger equation with the QCD Coulomb 
potential. To first approximation it describes a quark with momentum miv^ and an 
antiquark with momentum m2V^, where is the four- velocity of the meson. The spinor 
degrees of freedom for a non-relativistic pseudoscalar bound state are represented by 
the Dirac projector ^(1 -|-/)75. The non-relativistic approximation can also serve as a 
toy model for bound states of light (relativistic) quarks. We will in the following refer to 
"heavy mesons" as "5" (where we mean the realistic example of a Be meson, or the toy 
model for a Bq meson, with mi ^ m2) and "light mesons" as "n" (where the realistic 
example is rjc, and the toy-model application would be the pion, with mi ^ 1712). 



Definition of LCDAs for liglit pseudoscalar mesons 

Following [Hill we define the 2-particle LCDAs of a light pseudoscalar meson via 

(;r(/')|^i(3;)r^r5^2W|0) 
{7t{P)My)ir5q2ix)\0) 
{n{P)\qi{y)(7^yY5q2{x)\0) 

with two light-like vectors Zfi =yfi —Xfi and Pp, = P^ — m^/{2P -1)1^, with the usual 
gauge link factor [y, x] (Wilson line) understood implicitly. Here u = I —u denotes the 
light-cone momentum fraction of the quark in the pion, with ^;r(w) being the twist-2 
LCDA, while ^p{u) and ^crl") are of twist-3. For completeness, we have also quoted the 
twist-4 LCDA g7j:{u) which, like the 3-particle LCDAs, will not be considered further. 
All LCDAs are normalized to 1 , such that the prefactors in ([T]) are defined in the local 
limit X ^ y.ln the definition of ^^(a), we have included a factor 3/(Z) — 1), such that 
the relation between jdj^ and fiyt from the equations of motion (eom), 

flTt = I^K - {ini + 1^2) , (2) 
is maintained in D 7^ 4 dimensions. In the local limit the eom further imply 

f^n /"^ , X / \ 1 mi — m2 

1^71 = : , / duu(^p{u) = - + —- . (3) 

mi-fm2 Jo 2 IjAjt 

Notice that (I2I3I) hold for the bare parameters and distribution amplitudes. 

At tree level, and in leading order of the expansion in the non-relativistic velocities, 
the two quarks in the non-relativistic wave function simply share the momentum of the 
meson according to their masses, p^ ^ m,/ (mi -|-m2) P^- For "light mesons" this implies 

0;j(a) ~0p(m) ~g;r(M) =i 5(M-ao), (4) 
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where uq = m.\/{m\ +m2). Consequently, all positive and negative moments of the 
distribution amplitudes are simply given in terms of the corresponding power of uq. 
Notice that ~ at tree-level, and the corresponding LCDA ^^(m) can only be 
determined by considering the corresponding one-loop expressions (see 1121] )• 



Definition of LCDAs for lieavy pseudoscalar mesons 

We define the 2-particle LCDAs of a heavy pseudoscalar 5 meson following 161 [Till. 

aj8 
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where is the heavy meson's velocity, t = v ■ z and z = 0. Here /s is the 
(renormalization-scale dependent) decay constant in HQET. The Fourier-transformed 
expressions, which usually appear in factorization formulas, are given through 



0J=(O= / d(oe-"''^^{co), (6) 

J 

where 0) denotes the light-cone energy of the light quark in the B meson rest frame. 

Including a finite spectator quark mass m and the effect of the 3-particle LCDAs 
^A, as defined in [IS.] , the eom become 

— 2 

= (D-2) Tjr] r [^A(r?,^)-^v(r],^)], (7) 

Jo Jco-T] s dq 

which is trivially fulfilled at tree-level, where ^g{o)) ~ ^b{(o) ~ 5{(0 — m), and 
^y^(T], (^) = ff{(Xs). We have shown in [12] that this relation also holds after inclu ding 
Us corrections to the NR limit. (A second relation, which has been presented in [|l5l] 
and extended here to the case m 7^ 0, is found to be not valid beyond tree-level.) More- 
over, at tree level, the moments of the "heavy meson's" LCDAs can be related to matrix 
elements of local operators in HQET [6] . 



RELATIVISTIC CORRECTIONS AT ONE-LOOP 

The NR bound states are described by parton configurations with fixed momenta. Rel- 
ativistic gluon exchange as in Fig. [T] leads to modifications: First, there is a correction 
from matching QCD (or, in the case of heavy mesons, the corresponding low-energy 
effective theory HQET) on the NR theory, 

<^>M = 0i?^ + ^tf^ + ^(«.'). (8) 








FIGURE 1. Relativistic corrections to the light-cone distribution amplitudes. The dashed line 
indicates the Wilson line in the definition of the LCDAs. 



Secondly, there is the usual evolution under the change of the renormalization scale 
[16, t]]- In particular, the support region for the parton momenta is extended to < w < 1 
for light mesons and < (0 < oo for heavy mesons, respectively. 



Light mesons 

We first consider the leading-order relativistic corrections to the local matrix elements. 
We will focus on the case of equal quark masses (results for m\ ^ m2 can be found in 
[1121]). Our result for the decay constant. 
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is in agreement with fv^, and the result for and is consistent with the eom- 
constraint in (I2I3I) . using niji ~ mf -\-mf in the on-shell scheme. 

The remaining contributions to the NLO correction to the leading-twist LCDA contain 
an UV-divergent piece. 
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which involves the well-known Efremov-Radyushkin-Brodsky-Lepage (ERBL) evolu- 
tion kernel [1160, and a finite term. 
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Here the plus-distributions are defined as 
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TABLE 1. Convergence of (<^")7r moments (77 = as{p.) / as{m)). 



n 


2 


4 


6 


8 


10 


NLO {ji — m, in units of a^) 


0.333 


0.053 


0.019 


0.009 


0.005 


LL(T7 = l/5) 
LL(T7 = 1/25) 
asymptotic 


0.126 
0.173 
0.200 


0.048 
0.070 
0.086 


0.025 
0.038 
0.048 


0.015 
0.024 
0.030 


0.010 
0.016 
0.021 



An independent calculation of the leading-twist LCDAs for the rjc and 7/i// meson 
has been presented in [9]. Our result is not in complete agreement with their findings. In 
particular, we find that the LCDA quoted in [9] is not normalized to unity as it should 
be. 

At the non-relativistic scale, ji ~ m, the usual expansion of into Gegenbauer 

polynomials (the eigenfunctions of the leading-order ERBL evolution equations), does 
not converge very well, i.e. the Gegenbauer coefficients drop off slower than 1/n. A 
better characterization of the LCDAs at NLO is given in terms of the moments 

{^")^= f\u{2u-\Y^M, (13) 
Jo 

which are linear combinations of Gegenbauer coefficients of order < n . This corre- 
sponds to an expansion in terms of 5-function and its derivatives, 

Uu) = 2 £ {^")^^^ 5 W {lu - 1) . (14) 

Results for the first few moments {^")jt are shown in Table [IJ (Notice that the moments 
(^")ir. also receive corrections from sub-leading terms in the non-relativistic expansion, 
see (M-) 

Tabled] also contains the moments generated by ERBL evolution of the LO result dU 
for two values ofr} = (Xs{ll) / OCs{m) and in the asymptotic limit. To illustrate the change 
of the shape of under evolution, we employ a parametrization which is obtained 
from a slight modification of the strategy developed in UlSll . 

0;,(m)^— rdt{-\ntr-Un2u-\Jt'/')+f{2u-l-it'/')). (15) 
r[a,-lnftj Jo V / 

It involves three real parameters a> 0,b> 0,tc<l, which can be fitted to the first three 
moments {^^'^'^), and the generating function for the Gegenbauer polynomials, 

1 00 

Fig. [2ta) shows the evolution of the model LCDA as a function of u. For 7] = 1/5 the 
functional form still "remembers" the non-relativistic profile, while for 7] = 1/25 the 
LCDA gets close to the asymptotic form. 

The twist-3 LCDAs for the 2-particle Fock states are obtained in the same way as the 
twist-2 one. Details can be found in [fl2ll . In particular, all our results are in manifest 
agreement with the eom-constraint from ([3]). 




FIGURE 2. (a) Solid line: approximation for ^^{u) — 5{u~\ /2) in terms of the parametrization (fTSl l. 
and its evolution for T] = 1/5 (dashed line) and 7] = 1/25 (dotted line), (b) Evolution of the heavy- 
meson LCDA (jigico) = 5{co — m). The three curves (solid, long-dashed, short-dashed) correspond to 
77 = 1/2,1/5,1/10. 



Heavy mesons 

The calculation of the LCDA for a heavy meson goes along the same lines as for the 
light-meson case. However, important differences arise because the heavier quark is to 
be treated in HQET which modifies the divergence structure of the loop integrals. As 
a consequence, the evolution equations for the LCDA of heavy mesons [TJ differ from 
those of light mesons. 

Let us focus on the distribution amplitude ^g{(0) which enters the QCD factorization 
formulas for exclusive heavy-to-light decays. In the local limit we derive the corrections 
from soft gluon exchange to the decay constant in HQET, 
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Notice that the decay constant of a heavy meson exhibits the well-known scale depen- 
dence fl9']. The remaining NLO corrections to the distribution amplitude contain 
an UV-divergent piece, which involves the Lange-Neubert evolution kernel 7+((y, ft)',ju) 
[0]. As has been shown by Lee and Neubert [20], the solution of the evolution equa- 



tion can be written in closed form. The resulting evolution of 0^((y,/i) is illustrated in 
Fig. [2b), starting from 5{co — m) at the non-relativistic scale, for three different values 
of 7] = as{ll) / CCs{m) and taking (Xs{m) = 1. From the double-logarithmic plot we can 
read off the asymptotic behaviour of the LCDA for co ^ and co ^ oo. As argued on 
general grounds [|7|], ^g{(0) develops a linear behaviour for (O 0, whereas for co ^ oo 
the evolution generates a radiative tail which tends to fall off slower than 1 /(O at higher 
scales. 



The finite NLO correction to ^g{(0) reads 
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with an analogous definition of plus-distributions as in (fT2]) . In contrast to the light- 
meson case, the normalization of the heavy meson distribution amplitude is ill-defined. 
Imposing a hard cutoff Auv ^ and expanding to first order in m/A\jy, we derive 
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The last two expressions provide model-independent properties of the distribution am- 
plitude which agree with the general results in [20]. The two phenomenologically rele- 
vant moments in the factorization approach to heavy-to-light decays read 



Ob{{1) = Ab(jU) 
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A similar analysis can be performed for the LCDA 0g (ft)), for details we refer to [|12ll . 
In particular, we can read off the anomalous dimension. 
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which describes the evolution of 0^(co,ju) in the Wandzura-Wilczek approximation, 
where 3-particle LCDAs (and the light quark mass m) are neglected. Among others, 
7 is needed to show the factorization of correlation functions in SCET sum rules for 
heavy-to-light form factors [21]. Another new result are the first positive moments of 
^gico) as a function of the UV cutoff. 
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where the divergent pieces again are expected to be model-independent. 



SUMMARY 



Non-relativistic qq bound states have been used as a starting point to construct light-cone 
distribution amplitudes in QCD and in HQET. We considered relativistic gluon correc- 
tions at NLO in the strong coupling at the non-relativistic scale, as well as the leading 
logarithmic evolution towards higher scales needed in QCD factorization theorems. We 
also studied certain model-independent properties of light and heavy LCDAs, including 
new results for the LCDA ^'^{o},jl)m. HQET. 
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